Signatures of Fractional Quantum Hall States in Topological Insulators 
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The fractional quantum Hall (FQH) state is a topological state of matter resulting from the 
many-body effect of interacting electrons and is of vast interest in fundamental physics [J, 0] • The 
experimental observation of topological surface states (SSs) in three-dimensional bulk solids has 
allowed the study of a correlated chiral Dirac fermion system, which can host a single Dirac valley 
without spin degeneracy pj-Hol]. Recent theoretical studies suggest that the fractional quantum 
Hall effect (FQHE) might be observable in topological insulators [ll|, Q2]- However, due to the 
dominant bulk conduction it is difficult to probe the strong correlation effect in topological insulators 
from resistivity measurements [1, |||. Here we report the discovery of FQH states in Bi2Te3 from 
thermopower measurements. The surface thermopower is ten times greater than that of bulk, which 
makes possible the observation of fractional-filled Landau levels in SSs. Thermopower hence provides 
a powerful tool to investigate correlated Dirac fermions in topological insulators. Our observations 
demonstrate that Dirac topological SSs exhibit strongly correlated phases in a high magnetic field, 
and would enable studies of a variety of exotic fractional quantum Hall physics and other correlated 
phenomena in this newly discovered chiral Dirac system. 



Topological insulators (TIs) are a new class of quantum 
states of matter with topologically protected conducting 
SSs, arising from the topology of the bulk electronic band 
structure |L3Hd7j ]. There are two distinguishing features 
of topological surface states. One is the existence of an 
odd number of Dirac cones on each surface, and the other 
is the helical spin arrangement @, @, ■ Theoretically, 
the relativistic nature of Dirac fermions is believed to sig- 
nificantly modify the electron-electron interactions, with 
the possibility to produce more robust ground states at 
the n = 1 Landau level (LL) in TIs than in conventional 
two-dimensional electron systems [20H23T ] . The unique 
spin texture and the coexistence of non-insulating bulk 
states also raise the intriguing question of whether TIs 
may host exotic FQH states owing to the non-trivial 
Berry's phase [l(|, huge Zeeman energy (24|, and the 
screening effect from bulk carriers [Til ] . The potential re- 
alization of more stable non-Abclian FQH states in TIs is 
of practical interest for topological quantum computing 

mm. 

There have been magnctorcsistance measurements at 
n = and 1 LLs in (Bii-^Sba^Ses Q and at n = 4 and 
5 LLs in Bi2Sc2Te However, the sub-integer oscilla- 
tions at n = 1 LL in (Bii-^Sb^Ses can only be resolved 
in the second-derivative trace and their best linear fit in- 
tersects the filling- factor axis at instead of 1/2, incon- 
sistent with a Dirac spectrum. Such a property makes 
the exact origin of these oscillations unclear Q . On the 
other hand, the FQHE in TIs is theoretically precluded in 
the \n\ > 1 LLs ,11!], suggesting the features in Bi2Se2Te 
unlikely to be ascribed to the FQHE. So far, transport 
studies of FQHE in TIs have been limited to samples with 



mobility below 3,000 cm 2 /Vs and the SS conduction is 
susceptible to conducting bulk states. 

To explore the existence of FQH states in TIs, we 
present thermoelectric measurements on the Bi2Te3 crys- 
tals. The surface mobility of these crystals is up to 
14,000 cm 2 /Vs [![, comparable to the Hall mobility 
(30, 000 cm 2 /Vs) of high quality graphene where the 
FQHE has been recently discovered |26[. We first ex- 
amined the dependence of thermopower S xx on temper- 
ature T in both metallic and nonmctallic samples (Fig. 
1, insets). Though S xx shows a low-T peak in all these 
samples, the peak of the nonmctallic samples Ql and 
Q2 is significantly stronger than that of the metallic 
sample Ml. These observed peaks indicate the occur- 
rence of phonon-drag effect that is expected to appear 
at — 29 K in high purity Bi 2 Te 3 crystals [27|. It has 
been demonstrated that the phonon-drag thermopower 
from a two-dimensional (2D) conducting layer on a three- 
dimensional (3D) crystal can display giant quantum oscil- 
lations due to the phonon intra- and inter-LL scattering 
in the presence of a strong magnetic field. In such a 3D 
system, surface electrons are dragged by non-equilibrium 
3D phonons of the whole specimen, while in a purely 2D 
system such as graphene, electrons of a wavevector k can 
only interact with 2D phonons of a wavevector q < 2k. 
In addition, the bulk thermopower is considerably sup- 
pressed due to the existence of two types of bulk carri- 
ers with opposite signs Q. Therefore, we expect that 
the magneto-thermopower of SSs is orders of magnitude 
larger than that of bulk in the high field limit. The ther- 
mopower measurement thus provides a powerful tool to 
elucidate the nature of the topological SSs that is difficult 
to be probed by the conductance measurement. 
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FIG. 1: Magneto-thermopower in high magnetic fields, (a) 
The thermopower response — S xx versus 1/H in sample Ql 
at temperatures T = 3.5, 5, and 7 K. (b) — S xx versus 1/H in 
sample Q2. The insets show the T dependence of the ther- 
mopower profiles in samples Ql and Ml in (a) and in sample 
Q2 in (b). The arrows mark the sub-integer dips resolved in 
1/H > 0.06 T _1 . The dashed lines indicate the oscillating 
minima in the low field regime. 

Figure 1, a and b, show the thermopower response 
—S X x versus the inverse magnetic field 1/H in samples 
Ql and Q2, respectively. Large LL oscillations begin to 
emerge at H > 8 T and their amplitude becomes smaller 
as T decreases from 7 to 3.5 K. A prominent feature of 
these oscillations is that at 1/H < 0.06 T _1 sharp dips 
(black arrows) appear, with an aperiodic spacing smaller 
than the oscillating structure in the low field limit. 

It is illuminating to compare the oscillations in —S xx 
and the Shubnikov-de Hass (SdH) effect in the conduc- 
tance tensor G xx , which was confirmed to arise from the 
2D SSs in the previous study Q. As shown in Fig. 2, a 
and b, the extrema in —AS XX coincide with the extrema 
in AG XX . This occurs because both — S xx (S xx < for 
electron-like carriers) and a xx peak when the Fermi level 
(E-p) aligns with each LL, whereas vanish when E-p lies 
between LLs. Furthermore, we observe pronounced LL 
splitting near 1/H = 0.061, 0.102, and 0.142 T" 1 (gray 
dashed lines in Fig. 2b). This splitting indicates that the 
degeneracy is lifted between top (+) and bottom (— ) sur- 
faces. Similar effect has been seen in strained HgTe 3D 
TIs [lj|. Here, a weak Te composition gradient in Bi2Te3 
breaks the inversion symmetry and generates displaced 
Dirac points (see the inset of Fig. 2a). By cleaving the 
crystal into bulk samples with a thickness t — 20 ~ 100 
//m, we obtain slightly different surface carrier densities, 
which then leads to two sets of Landau filling factor v in 
one piece of sample (Fig. 2c). Hence, we can pinpoint 
the top and bottom surface index fields B„+ and B v ~ 
from the periodic spacing of strong (black dashed lines) 
and weak (grey dashed lines) minima in AG XX for sam- 
ple Q2. Similar results were observed in sample Ql, with 
its bottom surface cleaved adjacent to the top surface of 
sample Q2. As shown below, the extracted bottom sur- 



face carrier concentration in Ql is almost the same as 
the top surface carrier concentration in Q2. The Fermi 
energy levels relative to the Dirac point in samples Ql 
and Q2 are sketched in Fig. 2d. 

In terms of the carrier concentration n e on one surface, 
B v is related to v by 



where 4>q = h/e is the magnetic flux quanta, h Planck's 
constant, e the charge of electron, and 7 the filling factor- 
shift. A shift with 7 = corresponds to a conventional 
spectrum, whereas a deviation from the zero-shift with 
7 = 1/2 implies a Dirac spectrum. The 1/2 arises from 
the n = LL at the Dirac point. In the following, we la- 
bel the filling factors as v s , where s = ± indexes the top 
and bottom SSs. With the B~} identified in both — S xx 
and AG XX , we plot them against integers (triangles and 
circles in Fig. 3, a and b). The slopes of the linear- 
fit to the data yield the carrier concentration n e = 9.27 
(7.61) xlO 11 and 7.37 (5.92) xlO 11 cm" 2 , with the Fermi 
wavevector k F = 0.034 (0.031) and 0.030 (0.027) A~\ for 
the top (bottom) SSs in samples Ql and Q2, respectively. 
The linear-fit intercepts the v axis at 7 = 0.45 ± 0.02 in 
Ql and 0.67±0.02 in Q2, consistent with a Dirac disper- 
sion. Hence, we are again convinced that the 2D Dirac 
states give rise to the LL indexing shown in Fig. 2. Fur- 
thermore, the weak-field Hall anomaly provides an inde- 
pendent measurement of the average surface wavevector. 
The value of kp derived from surface Hall conductance 
is in reasonable agreement with the quantum oscillation 
analysis. 

We then illustrate the fine structure of — S xx in the 
range of 1 < v < 2 for samples Ql and Q2 (Fig. 3, 
c and d). We observed narrow, reproducible minima at 
the fields B v s (dashed lines) with v s = 5/3 ± . In ad- 
dition to the minima at v s — 5/3^, a valley- like struc- 
ture at v s = 9/5 ± is discernablc in Ql, and the dip at 
v s = 9/5 + becomes prominent in Q2 as the tempera- 
ture decreases. As seen in Fig. 3, a and b, the frac- 
tional Landau fillings lie on a straight line with the inte- 
ger ones. To examine the oscillating profiles of the most 
pronounced sub-integer structures, we plot the — S xx ver- 
sus the filling factor calculated as v = n^oB^ 1 +7 (Fig. 
4, a and b). The — S xx {y) traces obtained at various n e 
are almost overlapped, and their minima are all located 
around v = 1.67 ± 0.02 (traces are displaced for clar- 
ity). Both features strongly suggest that the observed 
high-field structures are associated with the FQH states 
at n = 1 LL. Based on the theory of thcomopower in 
the QHE regime, a minimum in a xx should accompany 
a minimum in thermopower. As shown in Fig. 2b (black 
arrows), two sub-integer dips are clearly resolved in Aa xx 
with v s close to ±5/3. Noting that for conventional 2D 
electron gas in the FQH regime, p xy should be quantized 
with {ve 2 /h)~ 1 at filling factor v. Unfortunately, the 
present of the bulk conduction channel does not allow us 
to measure exact Hall magnitude. Even in most resis- 
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FIG. 2: Quantum oscillations from integer Landau level, (a) 
The oscillating component of thermopower —AS X x, obtained 
by subtracting a smooth background from —S xx in sample 
Q2. The inset sketches the Dirac cone position on both sur- 
faces at the existence of a slightly Te composition gradient in 
bulk samples. The Fermi energy Ey is fixed across the bulk 
sample, (b) Longitudinal conductivity AG XX with a smooth 
background subtracted is plotted against 1/H at selected tem- 
peratures T = 0.7, 1.6, and 5 K in the same sample. Dashed 
vertical lines mark the minima in —AS XX and AG XX , which 
are consistent with integer v of top (+) and bottom ( — ) Fermi 
surface. The variation of Dirac point positions in two surfaces 
with a constant Fermi level Ep is shown as an inset in (a), 
(c) Sketch showing how the composition gradient induces the 
top and bottom Landau level displacement, (d) Illustration 
of the energy difference between the Fermi energy level Ef 
(short lines) and the Dirac point of top or bottom surface 
states in samples Ql and Q2. 

tive TIs, Hall quantization is obscured by dominate bulk 
contribution and no quantized Hall plateaus have so far 
been observed. 

In spite of the limited data set, we may roughly esti- 
mate the lower bound of the gap energy of the 5/3 state 
(A5/3) from the T dependence of the thermopower and 

resistivity minima, both of which scale as e _A / 2T . As 
shown in Fig. 2b and Fig. 3, c and d, the 5/3 + state per- 
sists until 5 K, indicating A 5 / 3 > 5 K. This value is more 
than an order of magnitude larger than the correspond- 
ing gap (A0/3) in the GaAs system with a much higher 
mobility [3C| . but comparable to the gap in graphene 
with a similar mobility [26|. This is not surprising be- 
cause the n = 1 LL in the topological SSs is a mixture 
of the n — and 1 LLs in non-relativistic systems. It 
makes the FQH states in the n = 1 LL in the Dirac sys- 
tem more robust than those in the GaAs system [ll|, US- 
UI, El, H3] ■ Compared with graphene, the surface states 
in Bi2Te3 has only one Dirac valley with no spin degener- 
acy, analogous to a completely four-fold degeneracy lifted 
graphene system where the FQH states do not mix be- 
tween spin- and valley-bands. Moreover, the helical spin 
texture of topological SSs and the presence of conductive 
bulk states may lead to enhanced effective Coulomb in- 
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FIG. 3: Identification of fractional quantum Hall states. The 
left panel shows the 1/H positions of the —S xx and AG XX 
minima versus the filling factor v in sample Ql (a) and sample 
Q2 (b). The black and red triangles represent the integer 
and fractional fillings, respectively. The right panel shows 
the fine structure of thermopower response —S xx versus 1/H 
for 1/H < 0.06 T" 1 in sample Ql (c) and sample Q2 (d). 
Identified minima correspond to the fractional-filling factors 
illustrated as the vertical dashed lines. Note that the minima 
at v — 5/3 + and 5/3" are much stronger than those at v — 
9/5+ and 9/5". 

teractions, rendering the FQH states even more robust 
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in TIs 

We next extract the surface thermopower S xx from 
the observed thermopower response. The measured ther- 
mopower tensor Sij can be expressed as the sum, 



k—x,y 



(2) 



where 



Pij is the total resistivity tensor, 

J Ql 
Dt, 



x,y ik kj 



with I — b or s, the bulk or surface ther- 
moelectric conductivity tensor, and the bulk or sur- 
face conductivity tensor. Since p xx 3> p yx , & h xx 3> cr%. y , 

for nonmetallic Bi2Te3 in the 
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high-field regime, S xx can be approximated as 



-p xx G\ 



St 



(3) 



where S xx is the bulk thermopower and G s xx the surface 
conductance. The bulk thermopower only gives rise to a 
featureless background. The p xx G s xx /t term can be ob- 
tained from the resistivity measurements. We find that 
the maximum magnitude of p xx G xx /t ~ 0.027 and 0.01 
in Ql and Q2, respectively. From Eq. ([3]), we can ex- 
tract the — S xx versus H in Ql and Q2 (Fig. 4, c and 
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FIG. 4: Characteristics of the v 
merical simulation of the surface thermopower. (a) 

Sxx as 

a function of the filling factor v = n e 4>oB~ +7 in sample 
Ql at 7 K. (b) — S X x versus v in sample Q2 at 7 K. In (a) 
and (b) curves are vertically offset for clarity, (c), (d) The 
surface thermopower — S% x obtained by removing a smooth 
background contributed from Sxx m sample Ql (c) and sam- 
ple Q2 (d) at 5 K. The green curve is the fit to — Sxx based 
on the Dirac fermion and 3D phonon interaction model. The 
H positions at the integer-filling factors are marked by dot- 
ted lines. The fractional structures resolved at v" < 2~ are 
indicated by arrows. 

d, red curves). The peak magnitude of —S% x is in the 
range of 0.5 — 2.0 mV K _1 , which is more than an order 
of magnitude higher than that of the bulk ~ 30 /iV K _1 
at 5 K. Unlike conventional 2D systems where the ther- 
mopower magnitude roughly displays a linear field depen- 
dence [28j . the surface thermopower at higher order LLs 
such as n = 4 is comparable or even greater than that 
of lower LLs (n = 3). This giant oscillating magnitude 



and the specific field profile of the surface thermopower 
can be understood within the scenario of the 2D Dirac 
electron and 3D phonon interaction. 

Because of the relativistic dispersion of topological sur- 
face states, the wave function fy n of a Dirac electron in 
the nth LL is the superposition of the nth and (n — l)th 
LL wave functions of a non-relativistic electron. The mix- 
ture nature of the wave function significantly modifies the 
electron-phonon matrix element in the n > 1 LLs, lead- 
ing to a thermopower profile different from an ordinary 
2D system. 

Using a general model given in rcf. 29 and the wave 
function for topological SSs [24j | , we numerically simulate 
the thermopower induced by the integer Landau quanti- 
zation in Ql and Q2, with the phonon mean free path 
treated as a fitting parameter (Fig. 4, c and d, green 
curves). We include the average LL broadening width 
F = 3.5 and 7 meV for Ql and Q2, respectively. The 
simple electron-phonon interaction model (Eq. S5) docs 
not capture the fractional features, as FQHE is not in- 
cluded in the model. However, it reproduces the index 
field position and the oscillation magnitude from the in- 
teger Landau quantization. This suggests that the ob- 
served giant integer Landau oscillations can be explained 
by 2D Dirac fermion and 3D phonon interaction. A more 
comprehensive FQHE framework is needed to model the 
fractional-filling states in thermopower response. 

By performing thermopower measurements, we have 
resolved fractional Landau quantization of SSs at v = 5/3 
and 9/5. The observed gap energy at the 5/3 state is ten 
times larger than that of the non-relativistic electron sys- 
tems. The demonstration of the FQH states in the topo- 
logical surface bands opens the door to future studies of 
fractional quantum Hall effect physics in the topological 
insulator, which is expected to display strong correlation 
effects between chiral Dirac fermions. 

The authors would like to thank N. P. Ong, F. D. Hal- 
danc, L. Fu, and C.-X. Liu for helpful discussion. 
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